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ABSTRACT 

Observations indicate that roughly 60% of the baryons may exist in a Warm-Hot In- 
tergalactic Medium (whim) at low redshifts. Following up on previous results showing 
that gas is released through galaxy mergers, we use a semi-analytic technique to es- 
timate the fraction of gas mass lost from haloes solely due to mergers. We find that 
up to ^ 25% of the gas in a halo can unbind over the course of galaxy assembly. 
This process does not act preferentially on smaller mass haloes; bigger haloes always 
release larger amounts of gas in a given volume of the Universe. However, if we include 
multi-phase gas accretion onto haloes, we find that only a few percent is unbound. We 
conclude that either non-gravitational processes may be in play to heat up the gas in 
the galaxies prior to unbinding by mergers or most of the baryons in the WHIM have 
never fallen into virialised dark matter haloes. We present a budget for stocking the 
WHIM compiled from recent work. 
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1 INTRODUCTION 

The baryon budget shows significant evolution from z ~ 
3, and results in a n apparent baryon defi cit today 
l|Fukugita et all 1 19981 : iFukugita fc PeeblesI |2004 ). At high 
reds hift, most of the baryonic mass is in the Ly -g for- 
est (iFukugita et al.lll998l : [Pukugita fc PeeblesI 120041 ) . while 
at low redshifts over half of the baryons are as yet unde- 
tected. The consensus is that the majority of the 'miss- 
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2006 ; iDave fc Oppcnhcimcr 200'^ with temperatures in the 



range 10 — 10 K - commonly referred to as the Warm-Hot 
Intergalactic Medium (WHIM). 

The immediate question is: how is this WHIM pro- 
duced ? Some form of mass and energy injection is essential 
to create this hot reservoir of gas; this form of feedback must 
both regulate the gas in galaxies and the metal content of 
the Intergalactic Medium (IGM). There has been much nu- 
merical work to incorporate various feedb ack mechanisms 
in an attempt to solve this pu zzle (e.g., 'Cen fc Ostrikeij 
Il999i : ,Nath fc Silbi200il ; iDave et al...2001. : ,Kang et al.,i2005l: 



ICen fc Ostrikeill2006l : iDave fc Oppenheimej[2007l . and ref- 
erences therein). Cosmological simulations seem to suggest 
that gravitational collapse during galaxy formation can pro- 
du ce and maintain the majority of the WHIM at 10^ — 10^ 
K (Cen fc Ostriker"l999"; 'Dave et al."l999'; 'Ccn fc O strike^ 
[2006; Croft ct al. 2001; Dave fc Oppcnhcimcr 200(1). Su- 
pernova feedback provides another avenue to generate the 
WHIM; for star bursts of 100 Mq per year, as much 
as 20% of the hot gas in a Milky Way mass galaxy 
can b e unbound jScannapieco et ahl l2006l : iKobavashi et al.] 
However, SNe feedback may be a self regulat- 
ing process, in that a starburst also heats the remain- 
ing gas and may damp the star formation rate, which 
in turn would reduce t he fraction of unbound gas (e.g. 
IScannapieco et"aL I l2008l ). Combining these effects, it is 
commonly thought that galaxies with host halo mass > 
10" Mq lose < 10% of their gas through SN feedback, 
while low mass ha loes may be enti r ely depleted of gas 
by this mechanism(|Yepes et al. 19971: Mac Low fc Ferraral 



ll999l : lEfstathioul200 0':'Scann apieco et a l."2006'). A third pos- 
sibility is that the radiation from an accreting supermassive 
black hole could p ower large-scale winds to blow mass out 
of th e galaxy (see Scannapieco fc OhI |2004|: Murray et al.l 
20051: iHopkins et all l2005al lbl l2006l : ICroton et all l2006l : 
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Siiacki et al.l |2007^ . Por a fixed amount of energy, all the 
non-gravitational feedback mechanisms are more effective in 
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low mass galaxies due to their shallower potential. However, 
observations suggest that low-mass galaxies are in general 
more gas-rich and are loss likely to have suffered a gas blow- 
out (iKannappan 2004; Gcha et al. 200^1 

In lSinha fc HoUev-BockelmannI l|2o'o9l . hereafter, SH09), 
we show that hot gas is driven into the IGM by galaxy merg- 
ers. The amount of hot halo gas lost depends strongly on 
the energy of the merger; it is possible for low mass galax- 
ies to retain their gas in this scenario during low-speed or 
distant encounters. However, SH09 only estimated the mass 
lost during a single merger. When all the mergers in the 
Universe are considered, this could heat and drive a signifi- 
cant portion of the total baryon budget into the WHIM. In 
principle, this process could join AGN and star formation 
feedback as a way to populate the WHIM, and we find that 
this method operates preferentially in more massive haloes. 
To estimate the total fraction of gas released by mergers, 
we construct a series of analytic halo merger trees using 
a publicly a vailablj^emi-analytic E xtended Press- Schechter 
(EPS) code l|Parkinson et al.ll2008h . 

In Section[2]we describe the theory of halo merger trees, 
in Section[3]we outline the experiments designed to track the 
gas ejected via galaxy mergers, in Section |4] we present the 
results for the halo gas ejected by this process and Section[5] 
contains the discussion. 



progenitors of mass Mi that make up a halo of mass K'h for 
a small step in redshift space of Az\ . This is written as: 



dN 



1 d/(Mi|M2)M2 



dzi (Ml < Ma). (2) 



dMi Ml dzi Ml 

Specifying a minimum mass resolution Mrcs allows us to 
compute the mean number of progenitors Np with mass 
Ml in a mass range A/rcs < M\ < A/2/2 via the following 
equation: 

fM2/2 .Ty 

The mass fraction F of the final halo AI2 that is accreted 
below Mres can be estimated from: 



F = 



dN Ml 



(4) 



A binary merger tree can then be constructed given 
A/2 and Z2. We used this technique to construct a set 
of twelve merger-trees, which we outline in Section [S] 
We have assumed a flat ACDM cosmology with Q,h ~ 
0.044, ^dni = 0.214, = 0.742, erg = 0.796 and /i = 0.719, 
consistent with the WMAP 5-year cosmology parameters 
(|Komatsu et al.ll2009l '). 



2 CONSTRUCTING MERGER-TREES 

Observations re veal that we l ive in l ow-density, A-dom i nated 
flat Univ erse (Ricss ct al.' 'l998'; 'Perlm utter et al.l 1 19971 : 
ISpergel"et al. .2007,; Komatsu et al. .2009|) . In such a Uni- 
verse, haloes form hierarchically, with smaller haloes form- 
ing early on and merging into larger structures at later 
times. This process of halo formation is dictated by grav- 
itational processes, and an analytic formalism yields the 
number density of haloes as a function of mass and redshift 
(|Press fc Schechteilll974l ). However, this does not constrain 
the merger rates for any given halo as a function of red- 
shift. To this end, the Press-Schechter formalism has been 
extended to calculate a m erger history of a ha l o in the form 
of a binary merg er tree (|Lacev fc Cold 1 19931 : lBoweJll99ll : 
iBond et al.lll99ih . These merger trees are computationally 
much less expensive than an N-body simulation, and are 
widely used to explore and constrain theories of galaxy evo- 
lution, black hole growth, etc. We use the technique here to 
estimate the gas lost to t he WHIM via galaxy me rgers. 

In the EPS model of [Parkinson et al.l (|2008l '). the con- 
ditional mass function /(M1IM2) gives the fraction of mass 
from a halo with mass A/2 at a redshift Z2 that was contained 
in a progenitor halo of mass Mi at a previous redshift zi : 



/(Mi|M2)dlnA/i 



2 af{5i^S2) 



exp 



1 (^1-^2)' 



2 (a? 



diner 
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d In Ml , 
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where Si and & represent linear overdensities for collapse 
at redshifts zi and Z2 and a = a{M). The derivative of this 
equation under the limit zi — > 22 yields the number A'' of 



3 METHOD 

As shown in SH09, the amount of gaflreleased by a galaxy 
merger depends on the mass ratio and the original gas con- 
tent of the haloes. To incorporate this effect within a merger 
tree we take the following approach: we seed each halo with 
a gas fraction (/seed) and assume a galaxy merger with a 
mass ratio greater than 77min unbinds a fraction of this gas 
(./unb). We also assume that as the halo grows by diffuse ac- 
cretion from the IGM, it also accretes gas at the universal gas 
fract ion, increasing the halo gas content. Recent simulations 
(see iKeres et all l2005l . l20od : lOekel et ahllioogl ) show that 
gas does not necessarily heat up to the halo virial tempera- 
tures; the majority of the haloes at low-z are only accreting 
cold gas. To estimate the effect of this multiphase accretion 
models, we divided the halo gas mass into hot and coldcom- 
ponents in accordance with Figure 3 of fKeres et al.l l|2009l ). 
After this partitioning, we follow the same procedure, ex- 
cept now we only unbind gas from the hot gas component. 
Table [1] outlines the parameters for the twelve experiments. 

We designed these twelve experiments to bracket the 
likely effect that galaxy mergers have on populating the 
WHIM. A reasonable upper limit is set by allowing even 
minor mergers (77min = 0.1) to unbind a flxed fraction 
(./unb = 0.1) of the progenitor gas mass (run Minorl). 
Our lower limit is set by seeding only the massive haloes 
(Mhaio ^ Mmin = 10^" Mq ) with gas at the universal gas 
fraction and allowing only major mergers (r; = 0.3) to re- 
lease gas (run MajorS). SH09 found that roughly equal-mass 
mergers can release up to 20% of their initial gas mass, and 
since the merger rate (per halo pe r redshift) is relatively flat 
from 0.3 < < 1.0 (see Fig. 8 in lFakhouri fc Mall200a ). we 
argue that rjmin ~ 0.3, funh = 0.1 is a good average scenario. 



^ |http : //star-ww.dur ■ ac .uk/- cole/merger _trees/| 



^ Since we do not model star formation in our semi-analytic ap- 
proach, we will use the terms gas and baryons interchangeably. 
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Table 1. The initial parameters for the twelve merger trees. 
Column 1 is the minimum merger ratio for gas to get ejected 
from haloes, column 3 is the fraction of the halo mass used to 
seed newly-appeared haloes (equal to fib/Cdm for all experiments 
other than MajorS and MinorS), column 4 is the fraction of the 
halo gas that unbinds during a merger, column 5 is the minimum 
halo mass that can retain gas and column 6 shows the entire mass 
range of haloes for which merger trees were made. 



^min 


Run 


fsccd 


funb 




Mass range 


[-] 


[-] 


[-] 


[%] 


[logMo] 


[logMe] 




Majorl 


0.21 


10.0 




8.0 - 13.0 




Major2 


0.21 


random 




8.0 - 13.0 


0.33 


Major3 


random 


10.0 




8.0 - 13.0 




Major4 


0.21 


10.0 




10.0 - 13.0 




Major5 


0.21 


10.0 


10.0 


10.0 - 13.0 




Major-Keres 


0.21 


10.0 




8.0 - 13.0 




Minor 1 


0.21 


10.0 




8.0 - 13.0 




Minor2 


0.21 


random 




8.0 - 13.0 


0.10 


Minor3 


random 


10.0 




8.0 - 13.0 




Minor4 


0.21 


10.0 




10.0 - 13.0 




Minors 


0.21 


10.0 


10.0 


10.0 - 13.0 




Minor-Keres 


0.21 


10.0 




8.0 - 13.0 



Haloes more massive than 10 Mq , representing groups or 
clusters of galaxies, can not be faithfully modelled using this 
binary galaxy merger mechanism and have been left out. 

The input parameters are the final halo mass, M2, the 
initial redshift, zi — 10 and the mass resolution, M^cs- We 
explore a range of final halo masses from M2 = 10* — 
10" Mq. We use 100 logarithmically spaced mass bins to 
create a merger tree for a specific M2 at the present epoch. 
To account for cosmic variance, we run 100 realisations of a 
fixed halo mass. Thus, overall we create 100 present day halo 
samples with 100 realisations for a fair sample of possible hi- 
erarchical merger histories of structure in the Universe. For 
each merger tree, we set Mros = M2 x 10"^. For lO" M© 
haloes, this value of Mres is comparable to the mass of an 
individual dark matter particle in our numerical simulations 
(SH09). We tested the effect of changing Mros, zi and the 
number of redshift levels and found that our choices produce 
convergent results for the estimation of the unbound gas. 

With these merger histories, we follow all mergers from 
z = 10 to z = that lead to a halo of mass M2, and eject 
a fraction of gas from the mergers with mass ratios greater 
than rjmin- The cumulative sum of the unbound gas produces 
the total gas released in assembling a particular halo. This 
yields the fractional gas lost by z = on a per halo ba- 
sis. We repeat this process for 100 realisations, which pro- 
vides the variance in the gas lost. We can find the total 
gas released in generating all haloes in the Universe by con- 
volving with the co-moving n umber density of those haloes 
at 2 = (|Warren et al.|[2006l ). Summing over the final halo 
masses yields the effect of halo assembly on populating the 
WHIM. 



4 RESULTS 

Figure [1] shows the redshift evolution of the cumulative gas 
mass lost from all haloes in a co-moving Mpc^ volume for the 
run Majorl. To generate Figure [1] we first take the mean 
of 100 realisations for the unbound gas mass in each red- 
shift step for each halo. This unbound gas mass is added up 
along the redshift track to yield the cumulative mass at each 
redshift step and then multiplied by the co-moving number 
density of that particular halo at z = 0. This is the cumu- 
lative co-moving density of the unbound gas for one halo 
mass. Repeating this process for the 100 final halo masses 
yields the individual tracks spanning the x-axis. Figure [T] 
shows that the most massive haloes unbind the most gas at 
all redshifts, in spite of their lower number densities. For 
example, the current number density in a co-moving Mpc^ 
of a 10^^ M0 halo is ~ 10** times smaller than for a 10* Mq 
halo; so the mass-density of the 10* Mq halo is an order of 
magnitude larger than the 10^^ M0 halo. This biasing to- 
wards higher mass is explained by the hierarchical assembly 
of haloes - more massive objects today undergo many more 
mergers in the pasl|f|. 

Figure [2] shows the redshift evolution of the unbound 
gas mass over the total baryon mass in all the haloes con- 
sidered in the merger tree. We find that 9% and 24% of the 
baryons can be ejected by mergers for the Majorl and Mi- 
norl runs respectively. The mass range of 10^" — 10^^ Mq 
and 10* - 10^^ Mq contain 39% and 52% of the total col- 
lapsed mass in the Universe respectively. Thus, the IGM 
pollution caused by the mergers presented in this paper can 
only refiect the history of at most half the total matter. If we 
assume that the same pattern holds true globally, then the 
fractions presented here (Figure[2ll can be interpreted as nor- 
malised by the total baryonic matter density of the Universe. 
Notice that the fraction of gas lost increases more rapidly 
with redshift for J^min = 0.1 - this is because 10:1 mergers 
occur more frequ ently than 3:1 (e.g.. iFakhouri fc Ma|[2008l : 
iGenel et al.ll2009l ). 

We can interpret Figure [2] in the following way: in the 
Majorl run, the convergence to 10% of the universal gas 
mass is tantamount to saying that the average halo under- 
goes one major merger in a Hubble time, since we set major 
mergers to release 10% of the gas mass. Likewise, the con- 
vergence of the Minorl run can be understood by noting 
that minor mergers (77 > 0.1) are 2 — 3 times m ore frequent 
than major merg ers (r j > 0.3, IFakhouri fc Mai (see bottom 
panel of Fig. 8 in l2008l )V Thus, the overall unbound density 
converges to ~ 20 — 30% for the Minorl run. 

In run Major-Keres with multiphase accretion, we find 
that only ~ 2% of the gas can be released due to mergers. 
Since the simulations of SH09 only included hot gas, we 
chose to unbind only from that phase. In the multiphase 
scenario, too much gas is in the cold phase and hence, can 
not be released during mergers. Even addi ng a mechanism 
to heat cold gas by major mergers (Eqn. 4 lCox et al.ll200j ) 
does not convert enough cold gas into a hot phase to be 
unbound later. If the haloes are only accreting cold gas and 
this gas can not be unbound from the haloes before heating 

^ There may be an additional effect from the higher gravitational 
potential en ergy (oc v'^^^^) involve d during mergers of massive 
galaxies (see [Johansson et al.|[2009h . 
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Figure 1. The cumulative gas mass lost from haloes by mergers in 
a co-moving Mpc^ volume as a function of halo mass and redshift 
for Majorl. Despite their larger number density, smaller haloes 
systematically lose less mass than the bigger haloes. The gas mass 
lost is obtained by taking a mean of 100 realisations. The pixels 
reflect the bin size in mass and redshift. This figure is summing 
the unbound gas mass along redshift and multiplying by the co- 
moving number density of the corresponding haloes at z = 0. 

it first, then the gas currently populating tiie WHIM may 
not have ever fallen into virialised iialoes. 

In a given merger tree, a fraction of unbound gas is 
released by mergers between small lialoes. To isolate the 
WHIM fraction (Table O Column 5)created during tiie as- 
sembly of only the massive galaxies, we run two sets of 
merger trees witii a lower mass limit of 10^° M0. Tabled 
shows that most of the unbound gas that is released comes 
during the formation of the massive galaxies. In particular, 
Major4, with only the massive haloes, produces nearly all of 
the unbound gas produced in the Majorl run. 

Although small haloes merging with massive haloes do 
not eject any gas, these minor accretion events increase the 
gas content of the remnant. This could potentially increase 
the amount of gas released by massive haloes in future merg- 
ers. However, if processes hke SN feedback evacuate the gas 
from low mass haloes, these low mass haloes can not increase 
the gas content of the massive haloes. To mimic this effect, 
we run two sets of merger trees with a lower mass limit of 
10^" M0 and only allow haloes larger than Afmin to contain 
gas. In this scenario (Major5 and Minor5), all small haloes 
are completely devoid of gas and therefore do not contribute 
to the gas mass of the big haloes. With this constraint, we 
find a WHIM fraction of ~ 3% and 8% for rjmin ~ 0.3 and 
0.1 respectively. 
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Majorl 

Majors 

Major-Keres 

Minorl 

\ 

\ 




Redshift 

Figure 2. The evolution with redshift of the gas mass lost as a 
fraction of baryon density in the Universe baryon content in the 
merger trees assuming universal gas fraction. If wc assume that 
the trend holds globally, then the Y-axis can be thought to be 
normalised by p^- The shaded region shows the l-rr deviation at 
each redshift from 100 realisations. 



Table 2. A census of the unbound gas at z = produced by 
galaxy mergers. Column 2 is the mass of unbound gas in a co- 
moving Mpc'' , column 3 is the ratio of unbound gas mass and the 
expected universal gas content (Mjiaio X ilb/f^dm)- Column 4 is 
the ratio of the mean unbound gas in a co-moving Mpc^ to the 
mean gas mass left the haloes. Column 5 shows the fraction of 
the WHIM generated; we assume that the WHIM contains 60% 
of all the baryons. The gas mass in a co-moving Mpc'^ volume for 
the haloes considered here is 7.6 X 10^ Mq. 



Run 




< Mg„„,„„i, > 


< M„„b.ga» > „ 

— ■ IWHIM 

< Mga,.ga, > 


[-] [10«Mo/Mpc3] [-] [-] [-] 


Majorl 


6.7 ± 0.4 


0.09 


0.10 


0.15 


Major2 


3.5 ± 0.2 


0.05 


0.05 


0.08 


Majors 


5.7 ± 0.3 


0.07 


0.09 


0.12 


Major4 


5.0 ± 0.3 


0.09 


0.10 


0.15 


Majors 


1.0 ± 0.1 


0.02 


0.04 


0.03 


Major-Keres 


1.7 ± 0.3 


0.02 


0.03 


0.03 


Minorl 


17.9 ± 0.4 


0.24 


0.32 


0.39 


Minor2 


9.7 ± 0.3 


0.13 


0.15 


0.21 


Minora 


15.2 ± 0.4 


0.20 


0.29 


0.33 


Minorl 


13.4 ± 0.4 


0.24 


0.32 


0.40 


Minors 


2.9 ± 0.1 


0.05 


0.13 


0.08 


Minor-Keres 


3.8 ± 0.3 


0.03 


0.06 


0.06 



5 DISCUSSION 

In this paper we show that a significant portion of the WHIM 
can be generated by gas ejected from galaxies during merg- 
ers. Our semi-analytic prescription shows that up to ~ 25% 
of the gas (assuming universal gas fraction) in haloes of mass 
10** -10" Mq can be ejected by mergers. Given an observed 
gas mass at 2: = 0, it is possible to infer the typical gas mass 
that was unbound from assembling that halo (column 4, Ta- 
ble [2]) For comparison with SN feedback, a quiescent Milky- 



Way type halo with star formation rat e of 1-10 Mpyr"^ 
would unbind 2% of the gas content (jScannapieco et aP 
I2OO8). We also find that multiphase gas accretion drastically 
reduces the amount of unbound gas from mergers, down to 
a few percent of the gas mass. In contrast with previous 
numerical work involving non-gravitational feedback, where 
the effects of mass loss are severe in smaller haloes, this 
merger mechanism unbinds gas preferentially from massive 
haloes. There is no selective unbinding of gas from dwarf 
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galaxies, in line with observational evidence suggesting that 
dwarf galaxies are more gas-rich and therefore may not hav e 
suffered a gas blow-out (Kannaopan 2004':l Geha et al 

This form of gravitational feedback may even play 
a larger ro l e in regulating the stellar mass function: 
iKeres et al.l ()2009t ) show that simulated galaxies exhibit 
a discrepancy w ith the observed stellar mass function 
l|Bell et al.l |2003| ) for both high and the low mass galax- 
ies. They conjecture that the key to solving this discrep- 
ancy is through a feedback mechanism that works for halos 
> 10^^ Mq - akin to our scenario. If the merger-ejection 
process is very efficient, then the current day haloes may be 
very gas-poor. It is conceivable that the current stellar mass 
and the gas fraction of a galaxy constrains the mean stellar 
mass and gas content of the galaxies of the past. 

Overall, we find that for our most reasonable scenario, ~ 
15% of the WHIM can be generated through galaxy mergers. 
If previ ous work on large-scale gravitational sh ocks proves 
correct l|Cen fc Ostrikej|l999l : [Pave et all 1200 ll ). - 66% of 
the WHIM can be attributed to gas that may have never 
fallen into a halo. In addition, recent studies have shown that 
roughly 20% can be produced via non-gravitational feed - 
back, such as SNe and AGN (e.g., ICen fc Ostri"keill2006l ). 
Therefore, with these three mechanisms to populate the 
WHIM, it may well be true that the baryon budget is bal- 
anced. 
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